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ABSTRACT 

Weakly interacting massive particles (WIMPs) are a viable candidate for the relic 
abundance of dark matter (DM) produced in the early universe. So far WIMPs have 
eluded direct detection through interactions with baryonic matter. Neutrino emission 
from accumulated WIMP annihilations in the solar core has been proposed as a signa- 
ture of DM, but has not yet been detected. These null results may be due to small scale 
DM density fluctuations in the halo with the density of our local region being lower 
than the average (~ 0.3 GeV/cm 3 ). However, the accumulated neutrino signal from 
WIMP annihilations in the Galactic stellar disk would be insensitive to local density 
variations. Inside the disk, dark matter can be captured by stars causing an enhanced 
annihilation rate and therefore a potentially higher neutrino flux than what would 
be observed from elsewhere in the halo. We estimate a neutrino flux from the WIMP 
annihilations in the stellar disk to be enhanced by more than an order of magnitude 
compared to the neutrino fluxes from the halo. We offer a conservative estimate for 
this enhanced flux, based on the WIMP-nucleon cross-sections obtained from direct- 
detection experiments by assuming a density of ~ 0.3 GeV/cm 3 for the local DM. We 
also compare the detcctability of these fluxes with a signal of diffuse high energy neu- 
trinos produced in the Milky Way by the interaction of cosmic rays (CRs) with the 
interstellar medium (ISM). These comparative signals should be observable by large 
neutrino detectors. 
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1 INTRODUCTION 

Determining the nature of dark matter and detecting it has 
remained an outstanding problem in astrophysics. For sev- 
eral decades, significant evidence in favour of a cold dark 
matter universe has accumulated from observations of galac- 
tic clusters and large-scale structure (Davis et al. 1985), su- 
pernovae (Perlmutter et al. 1999) and the cosmic microwave 
background (CMB) anisotropies (Bennett et al. 2003). It is 
believed that if WIMPs populate the halo of the Milky 
Way today, then they could be detected via direct WIMP- 
nucleon interactions. Experiments for direct detection are 
currently underway, e.g. PICASSO (Barnab-Heider et al. 
2005), CDSM (Akerib et al. 2006), SIMPLE (Alner et al. 
2005), NAIAD (Girard et al. 2005), and DAMA (Bernabei 
et al. 2008). Unfortunately, to date, no such direct interac- 
tions have been observed. However, these direct detection 
experiments have been used for placing stringent limits on 
the WIMP-nucleon cross-section. These limits are obtained 



assuming a value of ~ 0.3 GeV/cm 3 for the local DM density 
in the Milky Way halo (Barnab-Heider et al. 2005; Lewin & 
Smith 1996; Gates, Gyuk, & Turner 1995). 

WIMP-nucleon interactions in the Sun (Press & Spergel 
1985; Silk, Olive, & Srednicki 1985) or in the Earth (Freese 
1986; Gaisser, Steigman, & Tilav 1986; Krauss, Srednicki, 
& Wilczek 1986; Belotsky, Damour & Khlopov 2002) would 
cause the WIMPs to lose enough kinetic energy to become 
trapped in their respective potential wells, which could in 
turn create an observable neutrino flux from DM annihila- 
tion in the core. However, similar to direct detection exper- 
iments, no such neutrino signal has been observed. 

A crucial factor to consider in light of these null re- 
sults is the irregularities in the distribution of DM in the 
halo. Density fluctuations on sub-parsec scales allow for the 
possibility for our solar system to reside in a region of sig- 
nificantly lower density than the average (Taylor & Babul 
2005). If this is the case, then the actual WIMP-nucleon 
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cross-section could be higher than estimates which assume 
a local density equal to the average value of ~ 0.3 GeV/cm 3 . 
Furthermore, it could explain the lack of any direct detection 
of DM on Earth or any neutrino signal from the Sun. 

Extending this line of reasoning, we consider surveying 
the rest of the Galactic disk for an expected neutrino signal 
from WIMP annihilations in other stars. A large-scale sur- 
vey would alleviate any possibilities of lacking signals from 
DM voids due to regional density fluctuations. Estimates on 
high energy neutrino signals produced in various astrophys- 
ical processes has been a topic of interest for many years 
(for details, see Berezinsky & Zatsepin 1970; Berezinsky & 
Smirnov 1975; Stccker 1979). However, a new generation of 
neutrino telescopes is now opening a new window into the 
field. Large neutrino detectors (either built or planned) are 
capable of providing us with potentially ground-breaking 
observations in the fields of particle physics, astrophysics, 
and cosmology. Among the proposed detection techniques, 
under-ice and under-water optical Cerenkov detectors may 
prove to be particularly effective. This technique has already 
been implemented very successfully in the Lake Baikal ex- 
periment (Aynutdinov et al 2003), and in the AMANDA de- 
tector at the South Pole (Ackermann et al 2005). The much 
larger, km 3 -size detector IceCube (Ahrens et al. 2004) is 
currently under construction in the Antarctica, while other 
detectors (ANTARES, NEMO and NESTOR) are also being 
built in the Mediterranean Sea (Candia 2005). 

In this paper, as we extend the analyses done on the 
solar capture and annihilation of WIMPs to all the stars in 
the Milky Way disk, we show that for any given line of sight 
along the disk, the estimated neutrino flux is larger than 
the expected flux from WIMP annihilations elsewhere in the 
Galaxy or halo. Identifying point sources may be difficult to 
achieve, since they are likely to be too weak for a detector 
to receive an unambiguous directional signal. Hence, we fo- 
cus on detectible signals corresponding to unresolved diffuse 
fluxes. 

In order to do this properly, a sound knowledge of other 
possible diffuse neutrino signals is also important, such as 
those produced from cosmic rays. It is well known that at 
low energies, CR interactions in the atmosphere produce 
mesons, which decay into a neutrinos (Volkova 1980; Gaisser 
& Honda 2002). To date, this is the only type of neutrino 
signal identified so far. Another contribution to the high en- 
ergy diffuse neutrino fluxes is that resulting from the inter- 
actions of the CRs present all over the Galaxy with the ISM 
(Domokos, Elliott, & Kovesi-Domokos 1993; Berezinsky et 
al. 1993; Ingelman & Thunman 1996). The very low baryon 
density in the ISM (~ 1/cm 3 ) implies that essentially all 
mesons produced in such interactions decay without suffer- 
ing any attenuation (Candia 2005). Measuring and compar- 
ing these fluxes could potentially give valuable information 
about the distribution and behavior of DM and cosmic rays 
in the galaxy. 



2 DARK MATTER ANNIHILATION & 
NEUTRINO SOURCES 

Dark matter could exist in several forms. There have been 
several particle candidates proposed in the literature such as 
axions, neutralinos, and Kaluza-Klein particles, but speak- 
ing more generally here, we maintain the typical WIMP 
mass, M x , to be between a few GeV and a few TeV (for 
details see Jungman, Kamionkowski, & Griest 1996). In this 
work, we consider WIMPs in this mass range. If such a 
DM particle were in equilibrium with photons in the early 
universe, as the temperature decreased, a freeze-out would 
occur, leaving a thermal relic density. The temperature at 
which this occurs and the remaining DM density depends on 
the annihilation cross-section and mass of the WIMP. The 
cross-section is considered to be small enough to maintain a 
relatively slow and steady rate of dark matter annihilation 
even once large scale structure formation began and is still 
continuing within that structure today. 

The rate of neutrino production in WIMP annihila- 
tions is highly model-dependent as the annihilation fractions 
to various products can vary largely from model to model. 
We do not consider direct WIMP annihilation chanels into 
neutrinos, as they are considered rare or nonexistent, but 
many indirect channels exist. WIMPs typically annihilate 
into pairs of quarks or vector bosons. The quarks may then 
fragment into hadronic jets containing a large number of 
particles. The mass of the parent DM particle will effect 
what kinds of decay chains will follow from the annihila- 
tions. However, by setting the mass of the WIMP greater 
than the top quark, M x > M t , where M t = 175 GeV, calcu- 
lating the final neutrino outcome is greatly simplified. The 
reason is that the top quark lifetime is short, and almost 
100% of the time decays as t — > W. The W boson, like the 
t, decays before virtually any energy loss and creates neu- 
trinos with equal branching ratios (each about 10.5%) into 
ev e , fiv^, and tv t (Crotty 2002). Therefore, this is the only 
decay chain that will need to be considered. 

The two main sources of neutrinos from DM annihila- 
tion in the galaxy that we are considering here are from the 
general annihilations occurring in the halo, and for WIMPs 
that have been captured and have annihilated in the stars. 
Signals from general DM annihilation throughout the Galac- 
tic halo have been considered by (Finkbeiner 2005, Baltz & 
Wai 2004) , yet an enhanced signal ought to be present from 
DM annihilations inside the stars. In general, dark matter 
particles must be gathered in high concentrations to provide 
an observable flux of neutrinos. The potential wells of the 
Sun, Earth, stars or galactic centre are examples of regions 
where such concentrations may occur. Gravitational inter- 
actions bring a portion of the WIMPs in the Galactic halo 
and disk into the stars. For main sequence stars, we assume 
that the star is in hydrostatic equilibrium and is spherically 
symmetric, and that effects of rotation and magnetic fields 
are negligible (Lopes, Bertone, & Silk 2002). When WIMPs 
enter a star, they may interact with nuclei and lose enough 
kinetic energy to be trapped by the gravitational potential 
well and eventually annihilate with each other. 

There are two different channels by which WIMPs can 
scatter off nuclei in a given star: spin-dependent (axial) in- 
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teractions and spin-independent (scalar) interactions. The 
spin-dependent interactions refer to WIMP-on-proton inter- 
actions, while spin-independent interactions refer to WIMP 
interactions with the heavier elements in the star. Since an 
average star's content consists primarily of protons, the spin- 
dependent cross-section is typically several orders of mag- 
nitude larger than the spin-independent cross-section, and 
therefore, spin-dependent scattering dominates (Hooper & 
Silk 2004). 

As the WIMPs annihilate in the star's core, most of 
their decay products are subject to absorption and/or dif- 
fusion, but the neutrinos stand the best chance of escaping 
the star and propagating to earth relatively undisturbed. 
For a maximal neutrino signal, we assume at present that 
the WIMP capture and annihilation rates are in equilib- 
rium, Ta = Tc/2, (each annihilation destroys two WIMPs). 
Among the varieties of neutrinos produced, the muon neu- 
trinos, z/,,, are particularly useful for indirect detection of 
DM annihilation processes (hereafter 'neutrinos' specifically 
implies muon neutrinos). When the neutrinos collide with 
earth materials they produce muons, which have a relatively 
long range in a detector medium like ice or water. For more 
details on this detecting process of galactic neutrinos, see 
the results of (Abe et al. 2006; Kelly et al. 2005). 



3 CALCULATING THE NEUTRINO SPECTRA 
FROM WIMP ANNIHILATION IN STARS 

Assuming the Galaxy to be rotationally symmetric in the 
plane, the neutrino flux produced from WIMPs that have 
been captured by and annihilated in stars, $„„ , can be de- 
termined from the general expression 
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Here n cut <ri< ri max , where n represents the distance 
along a line of sight (LOS), ri cut is the cut-off distance from 
the Earth, which we take to be ~ 1000 AU and ri max is a 
great distance where the number density of stars in the disk 
is zero. We define r as the distance from the Galactic center 
to the source of the annihilation, geometrically related to n 
as 
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where I is the angle between the LOS to the Galactic center 
and the LOS to the source of the annihilation (n) for an 
observer on Earth. The emissivity at a point n is related 
to the dark matter density, pr>M, the stellar mass density, 
p*(r), and the emission rate per solar size star, L v , and is 
as follows, 
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In this expression p© = 6.5 x 10 25 g/cm 3 is the dark mat- 
ter density in our local region of the galaxy and /9dm (r) is 



the dark matter density at a point r, which we approxi- 
mate with a Naverro-Frenk- White profile (Navarro, Frenk, 
& White 1996) 
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where p — 5.6 x 10~ 25 g/cm 3 , and the scale radius, r B = 
20 kpc, (Finkbeiner 2005). The number density of stars in 
the midplane of the disk was estimated according to the 
following expression from (Efstathiou 2000) where 
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Here h d is the thickness of the galactic disk, 200 pc, 
and p,d = p e~ r ^ 3 ' 5kpc is the surface mass density of stars 
at a distance r from the galactic center. Considering lines of 
sights in the midplane, z = 0, we get 
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where p = 567Mq/pc 2 . In calculating the emission rate, 
L v , we account for all neutrinos within the energy range, 
-Emm < E v < M x , where E m i n is set by a typical detector 
threshold of 50 GeV (Crotty 2002). 
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where the emergent flux is 
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Here Ta is the annihilation rate of the WIMPs inside the 
star, ^ is the initial energy distribution of the neutrinos 
produced per annihilation, and P is the probability that a 
neutrino produced in the star's core will emerge from its edge 
without being absorbed by the steller material. Assuming 
that all the WIMPs annihilate at a rate equal to the rate at 
which they are captured, we define the annihilation rate as 
the product of the DM flux and the WIMP-nucleon cross- 
section: 



Ta = 4> x a* 

The DM flux, <j> x , is defined as (Barger 2002) 

1.2 X 10 7 _o _1 
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Assuming that the sun is an average size star, we may as- 
sume a reasonable capture cross-section for the sun cr© to 
be similar to the average capture cross-section for the stars, 
<r* , in the galactic disk. 
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where / = 10 is the focusing factor, i.e. the ratio of kinetic 
and potential energy of the WIMP near a given star, and 
1.2 x 10 57 is the number of nucleons in stars the size of our 
sun. The choice of parameters needed to compute WIMP- 
nucleon elastic scattering, a e \, is at best approximate and 
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are difficult to obtain with any generality as terrestrial ex- 
periments have indicated (Barnab-Heider et al. 2005). To 
this extent, we rely on the results of the various experi- 
ments, (Barnab-Heider et al. 2005, Gates, Gyuk, & Turner 
f 995, Lee et al. 2007) and take the most conservative cross- 
section estimates of WIMP-proton and WIMP-neutron in- 
teractions. We use a weighted-sum of the two, working out 
to ~ 0.5 — fpb which may vary slightly with WIMP mass. 
However, accounting for small scale DM density fluctua- 
tions, the WIMP-nucleon cross-section could be significantly 
higher than what we use here. 

The spectrum of neutrinos, is straightforward to 
calculate. Resulting from the decay chain t — > W — > u, 
we definine x as E/M x and integrate over the fragmenta- 
tion function of (Albuquerque, Hui, & Kolb 2001) to first 
compute the total number of hadrons, Nh, produced per 
annihilation. 
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where rj = Aqcd/M x , with the QCD scale set as Aqcd 
0.1 GeV, and the fragmentation function is 
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a = 15/16. The number of top quarks that are produced per 
annihilation is determined simply by the expression 
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From the number of top quarks produced per annihilation 
we can estimate the number of muon neutrinos that fol- 
low by assuming a branching fraction, B, of ~ 10% (Albu- 
querque, Hui, & Kolb 2001; Barger 2002). The number of 
top quarks and corresponding neutrinos produced per anni- 
hilation varies with the mass of the WIMP as indicated in 
Figure 1. 

The initial energy spectrum of the neutrinos produced 
per annihilation in the stellar core is: 
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100 is the normalization factor, B is the 



branching ratio, x — E/M x is the fraction of the initial 
WIMP energy given to the neutrino, and Mw = 80 GeV 
is the mass of the W boson. We also account for the possi- 
bility that when the tv t pair is produced by the W decay, 
the r also decays before losing much energy, producing a 
secondary around 18% of the time, (Crotty 2002) . 

The survival probability of these neutrinos must be con- 
sidered. Since a star is opaque to energetic neutrinos, the 
emergent energy spectrum per annihilation is not the same 
as the core spectrum. Some neutrinos which are produced in 
a stellar core suffer charged-current (CC) interactions and 
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Figure 1. (a) Number of top quarks (solid line) and neutrinos 
(dashed line) produced per WIMP annihilation as a function of 
WIMP mass, (b) The energy distribution of neutrinos produced 
per annihilation inside the stellar core (solid lines) and the final 
spectrum emitted from the stellar edge (dashed lines) for WIMP 
masses of (top to bottom) 500 GeV, 1000 GeV, and 5000 GeV, 
respectively. 



are rapidly thermalized as they make their way to the edge to 
escape. CC interactions remove a fraction of neutrinos above 
a certain transparency energy E^, where a(E] i ) J n(r)dr be- 
comes unity (Ritz & Seckel 1988). Thus, the probability of 
a neutrino escaping a star without interaction is given by 



(16) 



where E k = 150 GeV (Albuquerque, Hui, & Kolb 2001). A 
comparison of these spectra of neutrinos produced in the 
core and surviving to the edge of a star for different WIMP 
masses can be seen in figure 1. 
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4 NEUTRINO SPECTRA FROM DARK 
MATTER ANNIHILATION IN HALO 

Using the same general expression as equation (1) and in- 
tegrating over the same range, the neutrino flux from the 
'regular' DM annihilations in the disk and halo is calculated 
by setting 



where N — 3 x 10 , 7 = 1.63 is the spectral energy index, 
consistent with standard estimates for cosmic ray power- 
law spectra, (Berezinsky et al. 1993; Strong, Moskalenko, & 
Ptuskin 2007). We smooth over the inhomogeneities of the 
ISM density with a modest exponential function decreasing 
by a factor of ~ 2 from Galactic center to edge; ranging from 
~ 2 nucleons/cm 3 to ~ 1 nucleon/cm 3 , respectively. 
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where m x is the mass of a WIMP in grams and {av} — 
2 x 10~ 26 cm 3 /s is the DM annihilation cross-section, (as 
appears in Finkbeiner 2005). 



5 NEUTRINO SPECTRA FROM COSMIC RAY 
INTERACTIONS WITH ISM 

According to the standard picture, Galactic CRs are ac- 
celerated by the expanding shock wave fronts produced by 
the explosions of supernovae. This mechanism is capable 
of producing cosmic ray power spectrum with an index of 
~ 2.5 - 3 (Strong, Moskalenko, & Ptuskin 2007). As cos- 
mic rays propagate through the Galaxy, they interact with 
the ISM, which is a very low density, non-relativistic plasma 
constituted mainly by atomic and molecular hydrogen. The 
mesons produced in these proton-proton interactions decay 
before losing energy in secondary interactions. The muons 
produced from the meson decay then decay themselves and 
produce neutrinos. 

Given the uncertainties concerning the matter distri- 
bution in the Galaxy, one can only make a fairly simple 
model for it. The model used here is similar to previous 
ones (Domokos, Elliott, & Kovesi-Domokos 1993; Berezin- 
sky et al. 1993; Ingelman & Thunman 1996; Candia 2005), 
relatively standard, and should be adequate for our purpose 
of estimating the neutrino flux from CR-ISM interactions. 

In this short study of neutrino fluxes from CRs, inter- 
actions with stars are not taken into account. Even though 
CRs can indeed interact with the stellar material and pro- 
duce mesons, neutrinos and gamma-rays, these by-products 
must traverse the star such that the photon and neutrino 
fluxes would be significantly attenuated to a negligible level. 

In this scenario, in order to make a comparison with the 
neutrinos produced from WIMP annihilations we integrate 
over the same energy range here as well, _E m in < E < M x . 
The full expressions for calculating the neutrino flux can be 
found in (Candia 2005). Here we present the abbreviated 
expressions, i.e. the resulting emissivity in a parameterised 
form inspired by the analytic formalism of (Berezinsky et al. 
1993; Ingelman & Thunman 1996). Using the same general 
expression as equation (1) the neutrino flux from CR-ISM 
interactions in the disk is obtained with 
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where the values of 8 = 2.2 and 5 kpc were chosen to fit 
with the data provided by (Berezinsky et al. 1993). The in- 
teractions between CR particles and the ISM produce more 
p neutrinos than e~ neutrinos, however, by the time they 
arrive to Earth the neutrinos flavors are evenly distributed 
due to vacuum oscillations during propagation (Athar, Jez- 
abek, & Yasuda 2000). Since, the ISM and the CR densities 
are largest in the Galactic center region, there is a slight but 
noticeable anisotropy, being maximal in the direction to the 
Galactic center while being lower in the orthogonal planar 
direction as can be seen in figure 2 (thick solid lines). The 
comparison of the fluxes are presented in the next section. 



6 RESULTS 

In figure 2 the neutrino fluxes from WIMP annihilations and 
CR-ISM interactions are presented as a function of galactic 
longitude, I, for different WIMP masses and minimal en- 
ergy thresholds. Figure 2(a) compares neutrino fluxes from 
general DM annihilations, enhanced annihilations in stars, 
and from CR-ISM interactions for different values of M x . 
These spectra were all integrated from -E m i n = 50 GeV 
up to a maximum energy of M x . For each of the given 
WIMP masses, the neutrino signal from WIMP annihila- 
tions in stars is higher than general WIMP annihilations in 
the halo by at least an order of magnitude. For values of 
M x — 1000 GeV or lower, the neutrino signal from WIMP 
annihilation in stars is above the flux from CR-ISM inter- 
actions for ~ 5° or 10° from the galactic center. The flux 
produced from annihilations of heavier WIMPs would fall 
below the CR-ISM signal, primarily due to the significant 
attenuation the neutrinos suffer as they try to escape from 
the stars. 

Figure 2(b) compares neutrino fluxes from general DM 
annihilations, enhanced annihilations in stars, and from CR- 
ISM interactions for different minimal energy thresholds. 
These spectra were all integrated up to a maximum en- 
ergy of M x = 1000 GeV. Again, the neutrino signal from 
WIMP annihilations in stars remains consistently higher 
than general WIMP annihilations in the halo. For values 
of .Emm = 250 GeV or lower (well within detector con- 
straints), the neutrino signal from WIMP annihilations in 
stars is above the flux from CR-ISM interactions for ~ 5° 
or 10° from the galactic center. By -E m i n = 500 GeV the 
signal from CR-ISM dominates over the signal from WIMP 
annihilation. 

Fixing the viewing angle to be looking at the galactic 
center, Figure 3 compares neutrino fluxes from general DM 
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annihilations in the halo, enhanced annihilations in stars, 
and from CR-ISM interactions as a function of energy for 
different values of M x . Similar to Figure 2(a), these spectra 
were all integrated from i? m i n = 50 GeV up to a maximum 
energy of M x . In this figure, the neutrino signal from WIMP 
annihilations in stars is higher than general WIMP annihi- 
lations in the halo by at least a couple orders of magnitude. 
For designated values of M x = 500 GeV or M x = 1000 GeV, 
the neutrino signal from WIMP annihilation in stars is above 
the flux from CR-ISM interactions for energies up to several 
hundred GeV. These conditions of an enhanced WIMP sig- 
nal above the diffuse signal from CR-ISM interactions ought 
to be clearly observable within a fixed viewing angle and en- 
ergy range. 



7 DISCUSSION & CONCLUSIONS 

Detection of WIMP annihilation into neutrinos is a poten- 
tially promising indirect detection method. A small fraction 
of the neutrinos that reach the detector are converted into 
muons through charged current interactions (Gandhi et al. 
1996). These muons then propagate through the Cerenkov 
medium of the detector, where they are detected by pho- 
tomultiplier tubes distributed through the effective volume. 
For a review of high-energy neutrino astronomy, see (Halzen 
& Hooper 2002). As a muon propagates, it loses energy at 
a given rate and then may eventually fall below the energy 
threshold of the detector, typically 10-100 GeV for deep ice 
or water detectors. can be estimated as ~ E v /2. For ex- 
periments with energy thresholds like these, lighter WIMPs 
can be very difficult or impossible to detect for this rea- 
son. The neutrino-induced muon flux may be detected in a 
neutrino telescope by measuring the muons that come from 
stars in the direction of a line of sight in the Galactic disk. 
These neutrinos should have energies well above the energy 
of solar neutrinos and so should be easily distinguishable. 

At present, the AMANDA experiment is the largest 
high-energy neutrino telescope online, with an effective area 
of 50 000 m 2 and a muon energy threshold of about 30 GeV. 
It is currently taking data at the South Pole (Andres et 
al. 2001), and soon to be superseded by the much larger 
IceCube (Ahrens et al. 2004) , under construction there. Ice- 
Cube will have a somewhat higher energy threshold (50 - 100 
GeV) , but a full square kilometre of effective area. Also un- 
der construction, at higher latitudes in the Mediterranean, 
ANTARES (Blanc et al. 2003) will be significantly more sen- 
sitive to dark matter annihilations compared to AMANDA 
and IceCube due to its lower energy threshold (10 GeV) 
and larger effective area than AMANDA (100 000 m 2 ). And 
unlike the south pole telescopes, ANTARES will also be sen- 
sitive to neutrinos coming from the Galactic centre. These 
detector thresholds are well within the desginated WIMP 
masses chosen in this work. 

These kinds of experiments will have background events 
consisting of atmospheric neutrinos (Gaisser, Halzen, & 
Stanev 1995). In a given direction (up to the angular res- 
olution of a neutrino telescope), tens of events >100 GeV 
km~ 2 yr -1 are expected from the atmospheric neutrino flux. 
Fortunately, for detectors with very large- volumes and suffi- 




Figure 2. (a) Neutrino fluxes from CR-ISM interactions (thick 
solid line), general DM annihilations, and enhanced annihilations 
in stars as a function of galactic longitude for different values of 
M x . For M x = 500 GeV (dashed lines) and M x = 1000 GeV 
(dotted lines), the bottom and top curves represent the general 
and enhanced neutrino signals, respectively. These spectra were 
all integrated from -E m i n = 50 GeV up to a maximum energy of 
M x . (b) Similar neutrino fluxes as 2(a) but for different minimal 
energy thresholds. The thin solid and dashed lines represent the 
general and enhanced neutrino signals, respectively. These spectra 
were all integrated up to a maximum energy of M x = 1000 GeV. 
Read from top to bottom, all the lines represent the neutrino 
spectra with E m i n : 50 GeV, 250 GeV, and 500 GeV, respectively. 



cient statistics, it should be possible to significantly reduce 
this background or eliminate it entirely. 

Given the current limits of AMANDA, detecting a sig- 
nal from WIMP annihilations is not very strong. However, 
experiments with lower energy thresholds (ANTARES) and 
larger effective areas (IceCube) will greatly enhance the sen- 
sitivity and chances for detection in coming years. 

While this indirect detection method may appear 
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Figure 3. Looking toward the galactic center, neutrino fluxes 
from CR-ISM interactions (thick solid line), general DM annihi- 
lations, and enhanced annihilations in stars as a function of en- 
ergy for different values of M x . For M x = 500 GeV (dashed lines) 
and M x = 1000 GeV (dotted lines), the bottom and top curves 
represent the general and enhanced neutrino signals, respectively. 
These spectra were all integrated from E m i n = 50 GeV up to a 
maximum energy of M x . 



ing the neutrino survey up to higher energies. Given the 
high energies that cosmic rays can get to, they are capable 
of producing much higher energies than could be produced 
from DM annihilations, given the upper limits on the esti- 
mated WIMP mass. Also, at the relevant energy levels, the 
neutrinos produced from the CRs should be accompanied by 
a photon flux (Berezinsky et al. 1993), which would not be 
present from the DM annihilations occurring in stars. A dif- 
fuse gamma-ray flux above TeV energies from CR-ISM inter- 
actions could be measured with such ground-based detectors 
as CANGAROO, VERITAS, H.E.S.S. and Milagro, (Candia 
2005). In fact, the Milagro Gamma Ray Observatory has al- 
ready reported the first observations of a diffuse gamma-ray 
signal at TeV energies in the galactic plane (Atkins et al 
2005). Besides the pion decay channel, diffuse gamma-rays 
can also be produced by high energy electron bremstrahlung 
and inverse Compton scattering with the interstellar radia- 
tion field. For a detailed calculation of the photon flux pro- 
duced from CR-ISM interactions, see (Candia 2005). 

Being able to escape from very dense astrophysical en- 
vironments and travel vast distances undeflected and unat- 
tenuated, neutrinos offer precious information about DM as 
well as the high energy Universe. Such signatures will surely 
be subject to extensive experimental investigations in view 
of the new neutrino telescopes planned or under construc- 
tion. 



promising in some respects, it must also be kept in mind 
that the rate of neutrinos produced in WIMP annihila- 
tions is highly model-dependent. The sensitivity of a square- 
kilometre neutrino detector with a moderate muon energy 
threshold (50 GeV) will yield an observable signature in only 
a fraction of possible DM models, although such experiments 
are nevertheless an important probe. Accurately modeling 
the capture and annihilation rate of WIMPs by stars de- 
pends critically on DM distribution in the halo, and the 
WIMP-nucleon cross-section, which are all central to the 
prediction of the neutrino flux (and therefore muon flux) 
that is expected to be measured by neutrino telescopes. 

While most of our understanding of the properties of 
DM on subgalactic scales have come from numerical simula- 
tions of structure formation (De Lucia et al. 2004; Diemand, 
Moore, & Stadel 2004; Gao et al. 2004), they currently lack 
the resolution necessary to accurately estimate the level of 
substructure on the subparsec level. In which case, without 
such resolving power, it remains a definite possibility that 
the DM density in our local solar region is lower than the 
average DM density at the solar distance due to density 
fluctuations. We have thus chosen a very conservative value 
for the WIMP-nucleon cross-section based on estimates from 
terrestrial direct-detection experiments. Such a conservative 
cross-section would produce a rather modest capture and 
annihilation rate in stars. However, by surveying the stellar 
disk, we may see that the output flux from DM annihila- 
tions in distant stars could actually be higher than what is 
emitted from the Sun. 

The neutrino flux arising from CR-ISM interactions can 
be easily distinguished from the flux from DM by extend- 
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